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[l] Predicting the influences of dust on atmospheric com- 
position, climate, and human health requires accurate 
knowledge of dust emissions, but large uncertainties persist 
in quantifying mineral sources. This study presents a new 
method for combined use of satellite-measured radiances and 
inverse modeling to spatially constrain the amount and 
location of dust emissions. The technique is illustrated with a 
case study in May 2008; the dust emissions in Taklimakan 
and Gobi deserts are spatially optimized using the GEOS- 
Chem chemical transport model and its adjoint constrained 
by aerosol optical depth (AOD) that are derived over the 
downwind dark-surface region in China from MODIS 
(Moderate Resolution Imaging Spectroradiometer) reflec- 
tance with the aerosol single scattering properties consistent 
with GEOS-chem. The adjoint inverse modeling yields an 
overall 51% decrease in prior dust emissions estimated by 
GEOS-Chem over the Taklimakan-Gobi area, with more 
significant reductions south of the Gobi Desert. The model 
simulation with optimized dust emissions shows much better 
agreement with independent observations from MISR 
(Multi-angle Imaging Spectro Radiometer) AOD and MODIS 
Deep Blue AOD over the dust source region and surface 
PM 10 concentrations. The technique of this study can be 
applied to global multi-sensor remote sensing data for 
constraining dust emissions at various temporal and spatial 
scales, and hence improving the quantification of dust effects 
on climate, air quality, and human health. Citation: Wang, J., 
X. Xu, D. K. Henze, J. Zeng, Q. Ji, S.-C. Tsay, and J. Huang 
(2012), Top-down estimate of dust emissions through integration 
of MODIS and MISR aerosol retrievals with the GEOS-Chem 
adjoint model, Geophys. Res. Lett., 39, L08802, doi:10.1029/ 
2012GL051 136. 

1. Introduction 

[ 2 ] Mineral dust aerosol has important impacts on cli- 
mate change, air quality, and ocean biological systems 
[ Intergovernmental Panel on Climate Change, 2007]. How- 
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ever, significant uncertainties exist in mineral dust sources 
due to poor understanding of dust uplifting mechanisms and 
lack of in situ measurements over the desert region. Quanti- 
tative description of processes such as sandblasting and sal- 
tation bombardment in a chemistry transport model (CTM) 
requires knowledge of many parameters that can not be 
accurately characterized spatially and temporally, including 
surface wind speed, soil moisture, soil texture, surface cover 
type, and topography [Ginoux et al., 2001; Tegen et al., 
2002; Zender et al., 2003]. Not surprisingly, recent esti- 
mates in CTMs span from a few hundreds to over 4000 Tg for 
annual global dust emissions [Huneeus et al., 201 1] and vary 
by several orders of magnitude for emissions of dust events at 
regional scales [Uno et al., 2006]. An observation-based 
approach, therefore, is needed to reduce these large uncer- 
tainties in estimates of dust emissions. 

[ 3 ] In the last decade, in parallel with the advancement of 
in situ and remote sensing observation of dust aerosol, 
techniques have been developed to use these observations as 
constraints on dust sources. Cakmur et al. [2006] optimized 
dust emissions in a global model utilizing surface measure- 
ments of aerosol optical depth (AOD) and dust concentra- 
tions, providing an estimate of 1600-2500 Tg yr 1 . Dubovik 
et al. [2008] developed the GOCART adjoint modeling 
system for constraining global emissions of both fine and 
coarse mode dust from MODlS-based fine and coarse mode 
AOD, respectively. At a regional scale, Ku and Park [201 1] 
constrained soil emissions over East Asia by minimizing the 
difference between modeled and measured daily surface 
mass concentration of particulate matter less than 10 /itn in 
aerodynamic diameter (PM 10 ) in the downwind region 
through Bayesian least-squares analysis. Yumimoto et al. 
[2007] used a four-dimensional variational (4D-Var) sys- 
tem to optimize Asian dust emissions with constraints from 
ground-based lidar observations. 

[ 4 ] However, in the aforementioned studies, only a few 
independent observations far from the source region are used 
to validate estimates of dust emissions because ground- 
based observations of dust particles (particularly over the 
desert region) are sparse. Furthermore, while satellite AOD 
is a better constraint than surface measurements for esti- 
mating dust emissions as AOD represents the vertical load- 
ing of aerosols and thus the optimization of emission will be 
less sensitive to uncertainties in the simulated aerosol verti- 
cal profile, the discrepancy of aerosol optical properties 
between model and satellite can pose a challenge to attrib- 
uting the differences between satellite and modeled dust 
AOD to errors of dust emissions alone [Drury et al, 2008; 
Wang et al., 2010]. 
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Figure 1 . (a) True color image from MODIS/Aqua on 27 May 2008. The region within the rectangle of dashed line is the 
downwind plume region where the observational data are used. The region within the rectangle of dotted line is used for 
validation, (b) The prior mineral dust emission for the period of 25-28 May 2008. (c) Same as Figure lb but for optimized 
dust emission, (d) The ratio of the optimized to prior dust emission, (e) MODIS AOD at 0.67 /im aggregated to GEOS grid 
box and averaged for the period of 25-28 May 2008, which served as observational constraint for the adjoint inverse mod- 
eling of dust emission, (f) GEOS-chem total AOD at 0.67 //m simulated with prior dust emission averaged over satellite 
overpass time period (see text for details), (g) Same as Figure If but for the AOD simulated with optimized emission, (h) 
Ratio between Figures If and lg. 


[5] This study presents a new method to constrain dust 
emissions through integrated use of a chemical transport 
model (GEOS-Chem [Bey et al., 2001]) and its adjoint 
[Henze et al., 2007] with satellite-measured radiances. 
Although the adjoint optimization technique we use is sim- 
ilar to that of Dubovik et al. [2008] and Yumimoto et al. 
[2007], our study differs from the past studies in that: (a) 
the aerosol single scattering properties between GEOS- 
Chem and the satellite retrieval algorithms are consistent, 
and hence in essence it is the satellite-measured radiances 
that are used as constraints; (b) because of (a), the difference 
between satellite-based and model simulated AOD can, 
to first order, be attributed to model biases in the columnar 
loading of aerosols that is well suited to constrain the 
emissions; (c) optimization of dust emissions are conducted 
for every grid box as a function of time (e.g., on the daily 
basis); (d) the posterior dust emissions are evaluated spa- 
tially and temporally through comprehensive use of inde- 
pendent aerosol measurements from multiple satellite sensors 
and a field campaign. Specifically, the method is illustrated 
with a case study over East Asia for a dust event in May 2008 
(Figure la) during a joint China-U.S. dust field experiment 
[Huang et al., 2010]. The dust emissions in Taklimakan 
(hereafter TAK) and Gobi deserts are optimized daily using 
AOD retrieved from MODIS radiance [Wang et al., 2010] 
over the dark-surface region in China located downwind 
of the source region. In contrast, independent data sets are 
used for evaluating the optimization results over the dust 
source region, including MISR level-2 AOD [Kahn et al., 


2005], the MODIS level-2 bright-surface Deep Blue AOD 
[Hsu et al., 2006], sunphotometer-measured AOD, and 
ground-based measurements of PMi 0 from the field cam- 
paign during 2008 China-U.S. joint dust field experiments 
over Zhangye (39.08° N, 100.3° E, and denoted by green 
triangle in Figure la) [Huang et al., 2010]. 

2. GEOS-Chem Model and Its Constraints From 
MODIS 

[6] GEOS-Chem [Bey et al., 2001] is a global three- 
dimensional tropospheric chemical transport model driven 
by assimilated meteorological observations from the God- 
dard Earth Observing System (GEOS) of the NASA Global 
Modeling and Assimilation Office (GMAO). This model 
simulates the mass concentration of mineral dust in 4 size 
bins, together with sulfate, nitrate, ammonium, black and 
organic carbon in both hydrophilic and hydrophobic modes, 
and sea salt aerosols in both accumulation and coarse modes. 
The mineral dust entrainment and deposition (DEAD) 
scheme [Fairlie et al., 2007; Zender et al., 2003], combined 
with the GOCART topographic source function [Ginoux et 
al, 2001], is used in this study to calculate the first guess, 
or prior dust emission flux, as a function of wind friction 
velocity, soil moisture and soil particle size. 

[ 7 ] Dust emissions are constrained using MODIS AOD 
over the downwind dark-surface region (within the red 
dashed rectangle in Figure la) retrieved from an improved 
algorithm developed by Wang et al. [2010], Key to this 


2 of 6 


L08802 


WANG ET AL.: TOP-DOWN ESTIMATE OF DUST EMISSIONS 


L08802 


retrieval algorithm are: (a) a database of time-dependent 
surface 0.65 /zm and 2.1 fim reflectance ratios that are 
derived from samples of the MODIS dark-pixel reflectance 
data in low AOD conditions (e.g., “dynamic lower envelope” 
approach) for each GEOS-Chem model gridbox, and (b) an 
assumption that the simulated GEOS-Chem aerosol is unbi- 
ased in composition (but possibly largely biased in total 
mass), enabling use of the GEOS-Chem aerosol single scat- 
tering properties in the retrieval including the wavelength- 
dependence relationship needed for estimating AOD at 
2.1 //m from the AOD at 0.67 //m. With (a), (b), and a line- 
arized radiative transfer model that computes the top-of- 
atmosphere radiance and the sensitivity of radiance to the 
column AOD, Wang et al. [2010] is able to retrieve two 
unknowns (AOD at 0.67 /rm and surface reflectance at 
2.13 //m) from two known quantities (MODIS reflectances at 
0.67 and 2.1 qm). 

[8] Since the aerosol single scattering properties are con- 
sistent between the retrieval algorithm and GEOS-Chem, the 
retrieval essentially uses MODIS radiances (in the form of 
AOD) to scale the GEOS-Chem aerosol mass loading in 
each model grid column while keeping the aerosol compo- 
sition and aerosol vertical profile shape intact. Hence, this 
scale factor or the ratio of AOD between MODIS and 
GEOS-Chem (simulated with prior emission) is also the dust 
AOD scale factor that is needed for GEOS-chem simulated 
radiances to be consistent with MODIS radiances. Conse- 
quently, in cloud-free conditions, this (dust) AOD scale 
factor, twice per day respectively obtained from MODIS/ 
Terra and MODIS/Aqua retrieval, is used as a constraint for 
the GEOS-chem adjoint modeling. To ensure that the con- 
straint pertains to dust aerosols, only those (dust) AOD scale 
factors along the downwind dust transport path (e.g., 
retrieved dust AOD greater than 0.2) are used. 

[ 9 ] GEOS-chem simulated aerosol composition over Asia 
is shown by few studies to have low bias in black carbon 
(BC), and equivalent or larger underestimation of organic 
carbon (OC) mass and overestimation of sulfate aerosol 
mass [. Heald et al., 2005; Fu et al., 2012], which suggests 
that the mass fraction of highly absorbing (BC) and highly 
scattering (OC and sulfate) fine mode aerosols may have far 
less biases (as compared to the relative bias in OC mass 
only). Consequently, no significant biases are expected for 
(a) the GEOS-chem simulated fraction of coarse mode (dust) 
aerosol mass, and (b) the GEOS-chem simulated aerosol 
single scattering albedo. While (b) is important to ensure an 
unbiased retrieval of AOD, (a) supports that the GEOS-chem 
simulated dust AOD fraction is likely unbiased, both of 
which support the use of (dust) AOD scale factor derived 
from MODIS for constraining the dust emission. Admit- 
tedly, any model bias in modeled dust AOT fraction can lead 
to a corresponding bias (of the same sign) in the adjoint 
modeling results for dust emission. Further quantification of 
such bias requires a separate future study owing to the lack 
of aerosol composition data in China. 

3. Implementation of GEOS-Chem Dust Adjoint 
Modeling 

[ 10 ] The adjoint of the GEOS-Chem model was developed 
specifically for inverse modeling of aerosol (or their pre- 
cursors) and gas emissions \Henze et al., 2007]. The 
advantage of the adjoint approach is its ability to efficiently 


calculate model sensitivities with respect to large sets of 
model parameters, such as aerosol emissions at each grid 
box, which can serve as the gradients needed for inverse 
modeling of aerosol emissions. Here we implemented for the 
first time the adjoint of the GEOS-Chem dust simulation. 
Dust emissions are adjusted using a vector of control para- 
meters er that are the logarithm of emission scaling factors 
for dust emissions: a = ln(F/F a ), where F and F a are pos- 
terior and prior dust emissions vectors, respectively. The 
model response function, or the cost function, J, for the 
logarithmic dust emission scaling factor <r, is formulated as: 

J ( CT ) = “ c^rS^j/Ajo-) - c ofa ] 

z C ^n 

where c is the vector of model dust concentration (mixing 
ratio in the units of ppbv), H is the observation operator that 
transforms c to four-dimensional spatial and temporal 
observation space O, c obs is the vector of observed dust 
concentration described in Section 2, S Q b S is the observation 
error covariance matrix for dust concentration, 7 is a reg- 
ulization parameter, cr a is a vector of prior control para- 
meters, and S a is the error covariance matrix of er a . The 
adjoint optimization seeks to optimize cr by minimizing 
the cost function J through computing its gradient with 
respect to cr (i.e., V fT ./(rr)) iteratively in combination with 
the quasi-Newton L-BFGS-B algorithm [Byrd et al., 1995; 
Henze et al., 2007]. 

[ 11 ] Assumptions regarding the error in the observations 
and prior sources have important impacts on the optimiza- 
tion results. Based on assessments of the retrieved AOD 
(from MODIS) by Wang et al. [2010], we set the relative 
error for c obs as 20% and only consider model grid columns 
where dust mixing ratios are 0.2 ppbv and above because we 
are only interested in optimizing dust emissions. Similar to 
previous studies [Dubovik et al., 2008], S obs is assumed to 
have zero off-diagonal elements, and its diagonal values are 
the variance of each observation. To reflect the possibly 
large dust emission uncertainties in the model [Huneeus 
et al., 2011], S a is set to be diagonal with error of 200%. 
Lacking information to fully construct a physically repre- 
sentative S a , a regularization parameter (here 7 = 1000) is 
introduced to balance the contribution of model error and 
source error in the cost function, with a value selected using 
the “L-curve” technique [Hansen, 1998]. As a case study to 
illustrate the inversion of dust emissions, we first spin up the 
forward model simulation from 1 January 2008, and then run 
the forward and adjoint model iteratively for the period of 
20-31 May 2008 during which dust events occurred (25- 
28 May). The inverse modeling is conducted 5 days before 
the dust event to minimize the influence of initial conditions. 

4. Results 

[ 12 ] Figures lb— Id show the geographic distribution of 
the prior and posterior dust emissions, as well as the ratio of 
posterior to prior dust emissions during 25-28 May 2008. 
The posterior result in Figure lc is obtained when the cost 
function is sufficiently reduced (by ^60% after 8 iterations). 
As shown in Figure lb, the DEAD (or prior) dust mobili- 
zation scheme captures the overall distribution for the dust 
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Figure 2. (a) Scatter plot of GEOS-Chem AOD simulated with prior emission and the MODIS Deep Blue (DB) AOD over 
the dust source region for the period of 25-28 May 2008. The number of samples (N), bias, correlation (R), root mean square 
error (RMSE), and linear equation from reduced-major-axis regression are also shown, (b-d) Similar to Figure 2a but for the 
comparison of GEOS-Chem posterior (e.g., optimized) vs. MODIS Deep Blue AOD, GEOS-Chem prior vs. MISR AOD, 
and GEOS-Chem optimized vs. MISR AOD, respectively. In Figures 2a-2d, dotted and solid lines respectively show 
one-to-one fit and best linear fit between variables in x and y axis, while dashed line show the expected uncertainty in 
satellite-based AOD (i.e., ±20% for both MS1R MODIS Deep Blue), (e) time serial plot of GEOS-Chem simulated surface 
PM 10 concentration by prior (red line) and posterior (green) emission compared with the in situ measured PM 10 (black) 
over Zhangye station; also shown (on the upper right comer) are respectively the averaged value for these quantities over 
the same period, (f) Same as Figure 2e but for AOD time series with additional overlay of MODIS Deep Blue and MISR 
AODs (denoted respectively as blue and orange square). 


sources that are primarily located in the Takalaman desert 
and the Gobi deserts. While the prior and the posterior show 
consistent spatial patterns of dust emission over the Taka- 
laman region, the prior generally overestimates emissions by 
a factor of 2-3 over the Gobi desert. Overall, the optimiza- 
tion yields a 51% decrease in prior emissions, changing the 
total emission from 9.36 Tg to 4.55 Tg. The convergence 
between model and observation can be verified in that the 
model is in much better agreement with MODIS AOD 
(Figure le) after optimization in comparison with before 
optimization (Figure If). As expected, the adjustment of dust 


emissions in the optimization is consistent with the changes 
in GEOS-Chem AOD before and after optimization 
(Figure lh). 

[ 13 ] To validate the results of the inversion, we focus on 
the dust simulation over and near the source region where 
the dust has only been transported a short distance with 
minimal influence of precipitation and anthropogenic aero- 
sol. Figures 2a and 2d show the comparison of the GEOS- 
Chem AOD with independent measurements before and 
after optimization. The model AOD values in Figure 2 are 6- 
hour averages during 1:30-7:30 UTC that span the MODIS/ 
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Aqua and MISR/Terra overpass times, and the distribution 
of those AODs are shown in Figure SI in the auxiliary 
material. 1 When taking MODIS Deep Blue retrievals of 
AOD as reference, AOD values simulated from GEOS- 
Chem with posterior emissions have a better correlation 
(R = 0.88), smaller bias (0.01), and smaller root mean 
squared error (RMSE = 0.23) than the counterparts using 
prior emission (R = 0.77, bias = 0.16, RMSE = 0.39), which 
quantifies the improvements in the dust simulation with the 
optimized dust emission. Overall, the optimized GEOS- 
Chem AOD over 80% of grid boxes are within expected 
accuracy (100 ± 30%) of MODIS Deep Blue retrievals. 
Similar improvements are also found in the comparison of 
the model AOD to MISR AOD in Figures 2c and 2d, where 
bias and RMSE are reduced from 0.20 and 0.42 to 0.04 and 
0.15, respectively. 

[ 14 ] The improvement of the model dust emission is fur- 
ther reflected in the comparison of simulated near-surface 
dust concentration with ground-based observations. When 
compared to the ground measurements of PM 10 concentra- 
tion and AOD at Zhangye for the dusty time period, the 
modeled counterparts before the optimization of dust emis- 
sion correspondingly show a large positive bias or over- 
estimates by a factor of up to 5 (Figure 2e). This bias is 
significantly reduced in the simulation with posterior emis- 
sions. Moreover, the comparison of time series of AOD also 
shows that the modeled AOD after optimization has a much 
better agreement with MODIS Deep Blue and MISR retrie- 
vals (Figure 2f). 

5. Summary and Discussions 

[ 15 ] This study presents, for the first time, the integration 
of MODIS and MS1R data with the adjoint capacity for 
constraining dust emissions in the GEOS-Chem adjoint 
framework and conducts a case study that optimizes dust 
emissions in Taklimakan and Gobi deserts over East Asia for 
a dust event in May 2008. The model simulation with opti- 
mized dust emissions shows much better agreement with 
independent satellite and ground-based observations in the 
dust source region. This encouraging result confirms that 
using satellite observation combined with inverse modeling 
is an effective method for reducing the large uncertainties in 
estimates of dust emissions. 

[ 16 ] Indeed, estimates of dust radiative effect and forcing 
have been largely based upon modeling studies thus far 
[Forster et al., 2007]. Flowever, the semi-arid and arid areas, 
in global averages, are expanding with time due to global 
wanning, and hence, timely updates of dust source informa- 
tion, at finer spatial and temporal scales, is not only highly 
needed but also an outstanding challenge for chemistry 
transport models [. Denman et al., 2007]. Taking advantage 
of (i) adjoint modeling capabilities that can efficiently cal- 
culate the model spatial and temporal sensitivities of AOD 
with respect to model input parameters, and (ii) the high 
spatial and temporal coverage of satellite remote sensing 
of aerosols, the method developed in this study can be 
applied to over ten years of MISR and MODIS observations 
(with an estimated computational time of ^6 months if using 
a Linux cluster of 32 3.0 GHz CPUs), and hence has the 
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potential to improve and resolve the estimate of dust emis- 
sions at various spatial and temporal scales. Extension of the 
adjoint sensitivity tool developed here to examine the impact 
of specific mobilization parameters (such as the velocity 
threshold needed for dust mobilization) may advance the use 
satellite data to improve dust emission parameterization 
schemes. Consequently, such satellite-based constraints 
could improve our modeling of the dust radiative forcing, 
could potentially illuminate anthropogenic components of 
dust sources and loadings, currently estimated at 0-20% 
though values as large as 50% have been postulated [Forster 
et al, 2007], and could help refine understandings of longer- 
term feedbacks between land-use, climate and dust aerosol. 

[ 17 ] The results of this study are based upon analysis of a 
single dust event, and are therefore subject to model uncer- 
tainties regarding the relative bias in model simulated aero- 
sol composition and the aerosol size distribution, and to a 
lesser extent (because we use columnar AOD as a model 
constraint and select the days with little precipitation), the 
model errors in aerosol vertical profile and precipitation. 
While the use of model-simulated single scattering proper- 
ties in the retrieval of AOD from satellites strengthens the 
direct use of satellite radiances to constrain modeled AOD 
and aerosol emissions, satellite AOD products, for compu- 
tational expedience, may still be assimilated in regions where 
aerosol properties are dominated by one aerosol species (such 
as dust or smoke) and thus may be represented in stand-alone 
retrieval algorithms. However, in the long term, the use of 
satellite radiances to improve the aerosol model is the ulti- 
mate goal; this requires dedicated field experiments and in 
situ observations to provide evaluation of model improve- 
ments, as well as satellite missions that are dedicated to 
measurements containing information content beyond AOD 
(such as those with multi-angle multi-polarization to retrieve 
aerosol size and/or refractive indices) [Kahn, 2012]. 

[is] Nevertheless, as shown here with use of the satellite 
radiances, our finding of an overall 50% overestimation of 
dust emission in the model is consistent with Ku and Park 
[2011] (although there are some geographical differences 
due to the use of different versions of GEOS meteorology 
and thus different prior emissions), likely reflecting a sys- 
tematic bias in GEOS-Chem, and hence bearing important 
implications. A simple extrapolation of our results, using ID 
radiative transfer calculations with dust single scattering 
albedo of 0.96 in the mid-visible, implies a ~20% decrease 
of the dust direct radiative effect in annual and global aver- 
age. Further, if natural sources have been overestimated, then 
the anthropogenic fraction of the dust radiative effect (—0.5 
to +0.1 Wm“ 2 [Forster et al., 2007]) may have effectively 
been substantially underestimated. 
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